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ABSTRACT 

In this third paper in a series on stable magnetic equilibria in stars, I look at the stability 
of axisymmetric field configurations and in particular at the relative strengths of the toroidal 
and poloidal components. Both toroidal and poloidal fields are unstable on their own, and 
stability is achieved by adding the two together in some ratio. I use Tayler's (1973) stability 
conditions for toroidal fields and other analytic tools to predict the range of stable ratios and 
then check these predictions by running numerical simulations. It is found that while the 
poloidal field can account for no more than approximately 80% of the total energy, it can 
account for a very small fraction of the energy, i.e. that the toroidal field can be - and is likely 
to be - significantly stronger than the poloidal. Furthermore, the weaker the field, the weaker 
the poloidal component can be in relation to the toroidal. The implications of this result are 
discussed in various contexts such as the emission of gravitational waves by neutron stars, 
free precession, and a 'hidden' energy source for magnetars. 
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1 INTRODUCTION 

Magnetic fields are observed in various types of star which are con- 
sidered unlikely to harbour a suitable regenerative dynamo process 
because of the lack of convection. For instance, strong fields (300 
G to 30 kG) are observed via the Zeeman effect o n chemically- 
peculiar main-sequence A stars (the Ap stars; see iMathvsl 1200 ll 
for a review), as well as higher-mass O and B stars. These stars 
do contain a small convective core and a magnetic field produced 
inside it could in principle rise through the radiative envelope in 
the form of buoyant flux tubes, but the time-scale of this process 
is almost certainly too long for anything to be seen on th e sur- 
face within the star's lifetime (MacGregor & Cassinelli 2003]). The 
magnetic white dwarfs (WDs, with observed fields of 10 4 - 10 9 
G) have only weak surface convection, and neutron stars (NSs, 
with fields 10 8 - 10 15 G) no convection at all. All of these stars 
must therefore contain a stable 'fossil' magnetic field, either in- 
herited from the original molecular gas cloud or from the previ- 
ous stage_ofj^yohiti_onJe^ from the main sequence, 
see lWickramasinghe & Ferrarid J2005|). or neutron stars from a de- 
generate stellar core, see Ferrario & Wickramasing he (2006)) or 
left over from some kind of dynamo process at the time of for- 
mation, either during the pre-main-sequence phase, in the case 
of Ap stars, or during the convective proto-neutron star phase 
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dDuncan & Thompson! < fl992h ). To have survived since this time, 
a field must be stable on a dynamic (Alfven) timescal^]. 

It was suggested by Prendergast (1956) that a stellar mag- 
netic field in stable equilibrium must contain both poloidal (merid- 
ional) and toroidal (azimuthal) components, since both are unsta- 
ble on their own. It was then shown more rigorously by Tayler 
(1973) that any purely toroidal field configuration is unstable at 
at least some place in a star, and by Wright (1973) and Markey 
& Tayler (1973, 1974) that any purely poloidal field is unsta- 
ble. More recently, the properties of these instabilities of toroidal 
and poloidal fields have been looked at analytically and numeri- 
cally ( Sprui tll 19981, IBraithwaite & Spruifood Braithwaite 2006a, 
b, 2007. lBonanno &Urpinll2008h . 

Analytic methods have proven useful in demonstrating the in- 
stability of various equilibrium configurations, but have not been as 
useful in demonstrating the existence of any stable co nfiguration. 
Numerical simulations (Braithwaite & Nordlund 2006, hereafter 
Paper I; see also IBraithwaite & Spruitl l2004) showed that an arbi- 
trary initial magnetic field inside a non-convective star can evolve 
on an Alfven timescale into a stable configuration. A roughly ax- 
isymmetric configuration was found, of a mixed poloidal-toroidal 
twisted-torus shape as shown in fig. Q] Once formed it continues 

1 The solid crust of a NS can hold an otherwise unstable field in place, 
provided that the field strength is below some threshold ~ 10 13 G, but it is 
not clear whether the crust forms soon enough after the end of the proto-NS 
convective phase to prevent relaxation into stable equilibrium. 
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to evolve as a result of diffusive processes such as finite con- 
ductivity, on a much longer timescale. For instance, the diffusion 
timescale is ~ 10 10 years for an Ap star; in the case of a neu- 
tron star this timescale is much less certain and is a result of Hall 
drift and other processes as well as Ohmic dissipation. As the field 
evolves it moves outwards, passing quasi-statically through a se- 
ries of stable equilibria until, upon reaching the end of the series 
of axisymmetric equilibria, it changes to a non-axisymmetric equi- 
librium. These non-axisym metric equilibria are described in more 
detail in lBraithwaitd[2008l (hereafter Paper II), where it was also 
found that a non-axisymmetric equilibrium can be formed on an 
Alfven timescale directly, from somewhat different initial condi- 
tions. Essentially, the important difference is the central concentra- 
tion of the initial field and the fraction of flux connected through the 
stellar surface - a non-axisymmetric equilibrium can be formed di- 
rectly from an initial field whose energy is more 'spread out' rather 
than concentrated towards the centre of the star, and which has sig- 
nificant flux connection through the stellar surface. 



We are concerned here with only the axisymmetric class of 
equilibria. Since both toroidal and poloidal fields are unstable on 
their own, there is presumably some allowed range of ratios of 
the two respective field strengths; it is the principal purposes of 
this paper shed some light on what these stable ratios might be. 
The toroidal field is always confined to the interior of a star, since 
a toroidal field on or above the surface would require long-lived 
currents outside the star, so that we observe on the surface only 
the poloidal component. Therefore, it is either difficult or im- 
possible to get any direct observational constraint on the range 
of allowed ratios. There are however some interesting ways in 
which a toroidal magnetic field can manifest itself indirectly, which 
make the question of allowed poloidal/toroidal ratios a useful line 
of study. Firstly, in predominantly non-convective main-sequence 
stars (i.e. > 1.5Mq) it would be useful to know how much flux 
may be hidden below the surface, since this flux may be important 
during formation, eventually be visible on the WD, be responsi- 
ble for shaping the planetary nebula, affect the supernova explo- 
sion in some way, and affect the natal rotation periods of NSs and 
WDs via core-envelope coupling. In fact, this reminds us of an- 
other question: how much of the poloidal flux can be hidden be- 
low the surface? Certainly not all of it has to go through the stellar 
surface, but may be confined to the interior. This paper also sheds 
some light on this issue. In the context of the 'magnetars', highly 
magnetised NSs (dipole fi eld strength on surface 10 14 ~ 15 G; see 
Woo ds & Thompsorjliooi for a review), it would certainly be use- 
ful to know how much magnetic flux and energy could be 'hidden' 
in the interior of the star. These stars undergo soft-7-ray outbursts, 
releasing as much as 2 x 10 46 erg of magnetic energy in less than 
a second. A field of 3 x 10 14 G, if it fills the interior of the star, 
contains 2 x 10 46 erg, and since these stars appear to have a life- 
time of ~ 10 4 years and to undergo large outbursts perhaps once a 
century, it seems likely that the average field strength in the interior 
is significantly greater than that on the surface. 



Another way in which the poloidal/toroidal ratio may manifest 
itself is through its effect on the star's sha pe and moment of iner- 
tia. It has been know n for some time (e.g.lchandrasek har & Ferrm] 
ll953llWentzelll96lh that a poloidal field will make a star oblate 
and a toroidal field prolate, and obviously with a mixed poloidal- 
toroidal field it will depend on the ratio of the two. In general, such 



a deformed star should undergo torque-free pr ecessioifl; there is 
already some obs ervat ional evidence for this dCutler et alj I2003L 
IWassermanl 120031 and lAkgiin et al.1 120061) . If this precession is 
damped, then kinetic energy is minimised while conserving angular 
momentum and a prolate star will tend towards alignment of its ro- 
tation and magnetic axes, while in an oblate star the angle between 
the two axes will tend to 90°. In an Ap star, this damping pro- 
cess may or may n ot take on the order of a main- sequence lifetime 
dMesteletal1ll98lh but in a NS it may be much faster and a pre- 
dominantly toroidal field in a fast-rotating NS could lead to strong 
emission of gravitational radiation. These effects of the magnetic 
field on a star's moment of inertia will be looked at in more detail 
in a forthcoming paper (Braithwaite & Nissanke, in prep.) 

In section [2] the instability in toroidal fields is described and 
some predictions are made about the stability of mixed poloidal- 
toroidal fields, and the properties of instability in poloidal fields 
is reviewed. In section [3] analytic conditions are used to examine 
the stability of fields produced in simulations, and simulations are 
presented of the decay or otherwise of fields with various toroidal- 
poloidal ratios. I conclude, and discuss the results and their appli- 
cations in section|4] 



2 INSTABILITY IN AXISYMMETRIC FIELDS : 
ANALYTIC RESULTS 

In this section I review the nature and properties of instability in 
purely toroidal and purely toroidal magnetic fields in stars, as well 
as look at stabilisation of a predominantly toroidal field through the 
addition of a weak poloidal component and vice versa. 

2.1 The Tayler instability 

The st a bility of purely toroidal fields in stars w as examined by 
iTavlerl 119731) . who used the energy method of iBernstein et al.1 
( 1958) to derive necessary and sufficient stability conditions. Given 
a perturbation of the form 

in spherical polar coordinates (r, <f>, 0) he found that the stability 
conditions for m ^ 2 are less strict than for m — 1, so that if 
the goal is simply to distinguish between unstable and stable field 
configurations, these higher azimuthal wavenumbers need not be 
considered. The unstable modes can therefore be thought of as local 
in the meridional plane but global in the azimuthal direction. 

2.1.1 Necessary and sufficient stability conditions 

For stability against the m = mode we need the change in poten- 
tial energy, as given by the following integral, to be positive for an 
arbitrary displacement field £. In cylindrical coordinates (zu, <f>, z): 



SW oc / vjdvjdz{( f +a£l +b£ x €„ + c£i}, 



(2) 



where ( ) 2 is some function of £ which is positive definite, and 

where 



2 2 
P 9z 



dp 



(3) 



2 Purists may prefer the term 'nutation', although 'precession' occurs more 
frequently in the literature. 
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Tayler went on to show that it is impossible to satisfy all six of 
these conditions everywhere in the star, concluding that a stable 
equilibrium cannot be purely toroidal. 

A hand-waving explanation of the instability mechanism is 
as follows. One can imagine the toroidal field as a stack of field 
rings which exert pressure on one another. Magnetic tension ('hoop 
stress') and external pressure prevents the discs from simply ex- 
panding outwards and thus relieving this pressure, but it is energet- 
ically favourable for the rings to slip sideways out of the equilib- 
rium position, rather like the way an overloaded spinal column can 
result in a 'slipped disc'. The result is an m = 1 'kink' mode, as 
illustrated in fig. [2] 



Figure 1. The shape of the stable twisted-torus field in a star, viewed from 
different angles. The transparent surface represents the surface of the star; 
strong magnetic field is shown with yellow field lines, weak with black. 



dp dp _ 2pg z {p gza - 2Bl/zu) 2Ba dB^ 
dz dzu B, + ■yP zu dz 



Bl + jP 

dp {pg^-2Bl/mf _ 2Ba dB^ 1B% 
dzu B^ + jP zu dzu zu 2 

For 8 W to be positive, clearly it is sufficient (and can also be shown 
to be necessary) that the quadratic form is positive everywhere in 
the volume of integration. So the stability conditions are that 



a > , c > and b < 4ac 



(6) 



everywhere in the (zu, z) plane in the star. 

For stability against the m = 1 mode, we have a similar inte- 
grand consisting of one positive term and a quadratic where 

dp p 2 gl B% 2B 9B 



a = g-u. 



b = Ste-H- + 9- 



dzu 7P zu 2 zu dzu ' 
dp dp 2p 2 g^g z 2B$ dB^ 



dz dzu B\ + fP zu dz 



(7) 
(8) 



2.1.2 Unstable wavelengths and the stabilising effect of an 
added poloidal field 

The growth rate of the instability can be shown to be dTavleJl957h 



I'A 



Ba 



(10) 



In the case of an unstratified medium with infinite conductivity and 
zero kinetic viscosity, all vertical wavelengths (i.e. all n) are unsta- 
ble. However, finite conductivity damps the highest n modes faster 
than they can grow, resulting in a maximum unstable n. Conversely, 
stable stratification restricts the growth of the lowest wavenumbers. 
This is because, as we see from the continuity equation, the insta- 
bility results in vertical motions of order £ r ~ £hlr/lh, where £ r 
and £h are the displacements in the vertical (radial) and horizon- 
tal directions respectively, and l r and l\ are the length scales in 
the vertical and horizontal, so that l r = 1/n. The horizontal in- 
stability force per unit mass pushing a fluid element away from its 
equilibrium position is = <r 2 £h, and the vertical buoyancy force 
pushing the fluid element back to equilibrium is F r — N 2 £ r , where 
N is the buoyancy frequency. If the instability is to grow, we need 
therefore to have £h-Fh > £,rF r . Now, combining the se two limits 
we have the following range of unstable wavelengths dSpruitl 19991) 



a N 
> n > - — 



(ID 
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where r\ is the magnetic diffusivity. A toroidal field is stablised if 
the left hand side is less than the right hand side. These wa velength 
limits were confirmed numerically by Braithwaite] d2006al) . 

Adding a radial field component B r will also help to stabilise 
a toroidal field, as the horizontal displacements produced by the in- 
stability will bend the radial field lines, producing a restoring force. 
Using similar arguments to those in the previous paragraph, we can 
predict what field strength will be required for stabilisation. The 
(horizontal) restoring force per unit mass is ^B 2 /{lip) and, equat- 
ing this to the force from the instability, we have instability if 



£><s 



> n, 



(12) 



B r r sin 9B r 

using dlOb . This corresponds to a result in the case of an unstratified 

[ lasma from the unstable-modes and dispersion-relation method 
TavleJI 19571) . 

Taking numbers typical for a main-sequence star (in cgs units) 
r sin (9 ~ 10 u , p ~ 1, B T ~ B$ ~ 10 3 and r) ~ 10 2 , the term on 
the left of Jilt turns out to be 10 6 times greater than that in | |12I >. 
Similarly high ratios are found in WDs and NSs. This means that 
we can expect the stabilising effect of the poloidal field to dominate 
over that of magnetic diffusion. Therefore, ignoring magnetic diffu- 
sion and combining the effects of stable stratification and poloidal 
field, we have stability if 



N 



Bin 2 



<0, 



(13) 



which is true for all wavenumbers n if 

B% ^ 4B 2 pN 2 



7 4 < I 2 



(14) 



which can be rewritten in terms of the Alfven frequency uja 

B^jipOyfp) as 



N 2zu 



B,i 



(15) 



The second part of the left hand side will be roughly unity, so it is 
clear that the critical B r j B^ depends on the field strength in the 
star. In a main-sequence star, we have u\/N 2 ~ E/U, the ratio 
of mag netic to thermal ene rgies; in neutron stars the buoyancy is 
weaker ( Reisene ggej2008h . so that oj\/N 2 ~ 10 2 E/U. However, 
we expect E/U < 10~ 6 in a real star, so a relatively weak poloidal 
field should be sufficient for stabilisation in any kind of star. 



2.2 Instability of poloidal fields 

It is known that not only purely to roidal fields, but also 
purely poloidal fi elds are also unstable ( Markev & TayleiJ Il973l 
Braithwaite 2007). The instability begins in the region of the neu- 
tral line, the line where the poloidal component goes to zero. This 
instability is not unlike that in a toroidal field near the axis of sym- 
metry - the poloidal field increases in proportion to distance from 
the neutral line and it is the pressure in the direction parallel to 
this line, which the 'loops' are exerting on each other, which drives 
the instability. The direction of the stratification is the important 
difference between the two instabilities; displacements must be ap- 
proximately perpendicular to gravity, which in the poloidal case 
means that the loops move in a direction parallel to the star's axis 
of symmetry. This is illustrated in fig. [3] The growth rate of this 
instability is calculated in the same way as in the toroidal case, 
so it is given by the local Alfven frequency around the neutral 
line a) a = B v /{sy/p) where s is the distance from the neutral 




Figure 3. The form of the instability in a purely poloidal configuration. The 
left side shows the equilibrium and the right side the growth of a mode of 
particular azimuthal wavenumber. 



line. Adding a toroidal component (i.e. a component parallel to the 
neutral line) can stabilise the field, since the instability increases 
the length of the neutral line, thereby stretching the toroidal field 
line lying on it. The higher azimuthal wavenumbers are stabilised 
first, so that the m — 2 wavenumber is the first to become un- 
stable as the strength of the toroidal component is reduced. [Note 
that the m = and m = 1 modes are prevented by conserva- 
tion of momentum and angular momentum respectively.] In Paper 
II it was found that as an axisymmetric configuration diffuses out- 
wards towards the stellar surface, the toroidal component is lost 
into the atmosphere, weakening it relative to the poloidal compo- 
nent, and eventually the m = 2 mode appears and the configura- 
tion evolves quasi-statically along a non-axisymmetric sequence. 
The minimum-energy state of this non-axisymmetric equilibrium 
has comparable toroidal and poloidal components, which strongly 
suggests that an axisymmetric equilibrium cannot have a much 
stronger poloidal component than toroidal. A poloidal component 
which is stronger by some factor of order unity is possible, and 
the next section examines how large that factor can be. Note that 
this signals an asymmetry between the toroidal and poloidal cases, 
as there is no equivalent factor-of-order-unity argument giving an 
upper bound to the relative strength of the toroidal component. 

The strength of toroidal field required for stabilisation can be 
estimated with an analogous method to that described in section 
12.1.21 taking the wavelength to be 2nr n /m and arriving at the fol- 
lowing stability condition: 

B 



sB t 



m 
< —, 



(16) 



where s is the distance from the neutral line and r m is the radius of 
the neutral line. Note that B p is proportional to s in the neighbour- 
hood of the neutral line. Although this stability condition should be 
taken as rather approximate since curvature and other effects have 
been ignored, it shows again that stability against the (most unsta- 
ble) m = 2 mode requires the toroidal component to be at least 



comparable in strength to the poloidal component. IWrightl ( fl973h 
arrived at a similar stability condition (his eqn. 39): 

0.24S p (s ) < S t (IV) 

where so is some distance from the neutral line whose value is un- 
certain. Note that unlike in the case of strong toroidal fields in the 
previous section, the critical ratio here of toroidal to poloidal field 
strengths does not depend on the absolute magnetic energy in rela- 
tion to the thermal energy and buoyancy properties of the star. 

2.3 The effect of rotation 

In the previous sections the rotation of the star was not taken into 
consideration, but we may now look at its possible effect. First 
of all, we expect it to have an effect only in the 'fast' rotating 
regime, i.e. when the rotation speed is faster than the Alfven speed 
( Frie man & Rotenberg|l96(]|) . Applying this condition to real astro- 
physical objects, we see some fast rotators and some slow rotators. 
We find the ratio of the magnetic and rotation timescales to be 




where P is the rotation period, B is some average magnetic field in 
the interior of the star, and M and R are the stellar mass and radius. 
For 3M Q MS stars, O.8M WDs and 1.4M NSs, the numbers are 



— « 3800Prod55kG OT 18-Pid 1 ^ 1 or OmP^B^ (19) 

U)A 

where the periods are in units of 10 days, 1 day, 10 seconds and the 
field strengths 5kG, 10MG and 10 15 G - typical values for mag- 
netic MS stars, WDs and NSs (magnetars) respectively. Evidently 
the typical Ap star is in the fast rotating regime; however there are 
a few known examples with comparable field strengths but rotation 
periods of 30 — 100 y ears and above , just putting them into the 
slowly rotating regime (Mathys 2008). There is also much varia- 
tion amongst the magnetic WDs: a group exists with P ~ lOOyr 
and B ~ 300 MG while some others have P ~ few hours and 
B ~ IMG (Ferr ario & Wickrama singhe 2005f), so this class of star 
contains both fast and (very) slow rotators. The magnetars show 
much less variation and all are slow rotators; the radio pulsars on 
the other hand are fast rotators. 

In the fast rotating stars, the rotation will certainly have some 
effect on the growth rate of any instability. In the rotating frame, the 
Coriolis force acts perpendicularly to any unstable displacement, 
giving rise to epicycles; in the absence of any diffusive mechanisms 
this causes the fluid elements to come back to their original posi- 
tions and stabilisation results. This has been seen in simulations 
of the instability in a toroidal field dBraithw aite 2006a). Diffusion 
though damps the epicycles and the instabilit y returns, albeit with 
the growth r ate reduced by some factor n/u>A dPitts & Tavled 19861 
Spruit 19991) . In the case of a st rong poloid al field, a similar effect 
has been seen in simulations (Braithwaite 2007). It seems likely 
then that rotation does not affect the upper and lower bounds on 
E p /E. 



3 STABILITY OF AXISYMMETRIC FIELD 

CONFIGURATIONS: NUMERICAL METHODS 

Having derived some stability criteria for a mixed poloidal-toroidal 
field, it remains to see what this means for the stability of a global 



Axisymmetric magnetic fields in stars 5 

magnetic field configuration. Clearly, a globally stable configura- 
tion needs to be locally stable against the Tayler instability at each 
point in the star, as well as being stable in the neighbourhood of 
the neutral line. There are two obvious ways to proceed. The first 
is to find a global configuration and then check that a local analysis 
predicts stability at every location in the star for both types of in- 
stability. The second way is to construct a global configuration and 
then numerically follow its evolution in time. 



3.1 Constructing an axisymmetric field 

The basis for constructing axisymmetric field configurations whose 
stability we can probe will be the result of numerical simulations 
similar to those performed in Paper I, to where the reader is referred 
for a fuller account of the setup of the model; a brief outline is given 
here. 

The code used is the STAGGER CODE dNordlund & Galsgaardl 
1 19951 iGudiksen & Nordlund 120051) ■ a high-order finite-difference 
Cartesian MHD code which uses a 'hyper-diffusion' scheme, a 
system whereby diffusivities are scaled with the length scales 
present so that badly resolved structure near the Nyquist spatial 
frequency is damped whilst preserving well-resolved structure on 
longer length scales. This, and the high-order spatial interpolation 
and derivatives (sixth order) and time-stepping (third order) in- 
crease efficiency by giving a low effective diffusivity at modest res- 
olution (128 3 here). The code includes Ohmic and well as thermal 
and kinetic diffusion. Using Cartesian coordinates avoids problems 
with singularities and simplifies the boundary conditions: periodic 
boundaries are used here. 

The simulations model the star as a self-gravitating ball of 
ideal gas (7 = 5/3) of radius R in hydrostatic equilibrium with 
radial density and pressure profiles obeying the polytrope relation 
P oc p 1+ ~ , with the index n set to 3 here, since this gives stable 
stratification and is a reasonable approximation to an upper-main- 
sequence star. It seems unlikely that a different EOS, for instance 
that of a neutron star, will make even much quantitative difference 
to the results. The important point is the stable stratification - the 
issue of magnetic equilibria in a non-stably-stratified star will be 
explored in a forthcoming paper. 

A boundary between the star and the surrounding volume is 
produced in the following way. In reality, we expect the field out- 
side the star to be not only force-free (as is the usual approxima- 
tion when looking at magnetic loops in the solar corona) but also 
curl-free, since a very tenuous medium will not sustain long-lived 
currents. Numerically, a potential (i.e. current-free) field is tricky to 
produce, but fortunately the same effect can be produced by adding 
a high magnetic diffusivity to the volume outside the star, and in 
practice the diffusivity chosen is the maximum allowed without 
having to reduce the timestep. This causes the field outside the star 
to relax fairly rapidly to a curl-free configuration. Also, the gas out- 
side the star is hot, increasing the scale height and thus stopping the 
density from falling to greatly towards the edges of the computa- 
tional box, preventing high Alfven speeds and the lower timestep 
they would cause. The star is given an initially random magnetic 
field containing energy at all length scales down to a limit of a few 
grid-spacings, and the MHD equations are integrated in time to fol- 
low the evolution of the field. Within a few Alfven crossing-times, 
a stable equilibrium is reached. In the case where the field is less 
centrally concentrated, the equilibrium is non-axisymmetric, con- 
sisting of twisted flux tubes meandering at roughly constant depth 
under the surface of the star. This case was examined in detail in 
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Paper II. In the case of more centrally-concentrated (i.e. deeply- 
buried) initial fields, an approximately axisymmetric field forms. 
Ignoring the small deviations from axisymmetry, there appear to be 
two basic degrees of freedorrjf]: the concentration of the field into 
the centre of the star, which can be parametrized as r n , the distance 
between the axis of symmetry and the neutral line, and the poloidal 
fraction of total energy E p /E. In this two-dimensional parameter 
space lies an area of stability through which the star slowly moves 
in time, as a result of diffusive processes such as finite conductivity, 
to ever increasing r n and the eventual end of its stable axisymmet- 
ric life, as described above and in Paper II. It is the aim here to find 
the boundaries of that area of stability. 

The first step is to produce a stable field in a simulation (i.e. 
run with arbitrary initial conditions for a few Alfven crossing times 
until the field has settled down into an equilibrium) and axisym- 
metrise it, using an axis defined by J r x HdV. Although the fields 
produced in the simulations are already approximately axisymmet- 
ric, perfect axisymmetry simplifies the stability analysis. Also, the 
small asymmetry between the two hemispheres is removed so that 
the field is symmetrical about the 2 = plane. Now, this sym- 
metrised star can be put back into the simulation and evolved in 
time (the 'fiducial simulation'), where it slowly diffuses outwards, 
the value of r n rising as it does so. Its stability can then be ex- 
amined at various points along this diffusive evolution track, both 
analytically, using the methods described in Sect. [2] and numeri- 
cally, by changing the relative strengths of the poloidal and toroidal 
parts and using that as the initial conditions for a new simulation. 
In this way, the boundaries of the stable area in r n - E p /E pa- 
rameter space can be found. [In the simulations described here, the 
ratio of magnetic to thermal energies in the star E/U = 1/400 
unless stated otherwise.] In fig. [4] we see the axisymmetric field at 
the times t/r A = 5.2, 69.7 and 97.0 when r n /R = 0.33, 0.47 and 
0.58. Time is measured in units of the Alfven crossing time, defined 
here as ta = Ry/ M/2E where E is the total magnetic energy and 
M is the mass of the star. Note that it can be seen in the figure that 
the contours of vjB^ are parallel to the poloidal field lines. This 
result can be derived analytically jMestell 1 1 96 ll iRoxburghl 1 1 966h 
from V • B = and from the recognition that the azimuthal com- 
ponent of the Lorentz force must be zero everywhere because in an 
axisymmetric equilibrium it cannot be balanced by gravity pg or 
the pressure gradient —VP. 

In fig. [5] are plotted various quantities from the fiducial run 
which change in time as the magnetic field evolves diffusively: the 
neutral line radius r n , the poloidal energy fraction E p /E and the 
fraction of the poloidal flux which breaches the surface of the star 

$surf/$p. 

First, the case where the toroidal field dominates will be exam- 
ined, and then in section [33l the poloidal-dominated case is looked 
at, in order to find both boundaries of stability, i.e. both upper and 
lower bounds to E p /E. 

3.2 Stability of a predominantly toroidal field 

In this section I examine the stability of fields with low E p /E, 
firstly with a combination of Tayler's stability conditions and equa- 
tion d!4t . and secondly with simulations using different values of 
Ep/E as initial conditions. 

3 There are other, more subtle degrees of freedom concerning the exact 
shape of the field lines, but they appear to be less important, and are beyond 
the scope of this study. 




Figure 4. Projection onto the meridional plane of stable magnetic field con- 
figurations with r n /R = 0.33, 0.47 and 0.58. The yellow semicircle is 
the surface of the star, and the black lines and blue shading represent the 
poloidal and toroidal components of the field, respectively. The contours of 
the toroidal pail are actually contours of ■ujB^; it is clear that this quantity 
is roughly constant on poloidal field lines. The poloidal lines plotted are 
separated by equal quantities of poloidal flux. 
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Figure 5. The neutral line radius r n , the poloidal energy fraction E p /E 
and the fraction of poloidal flux which breaches the surface of the star 
'^'surf/'I'p against time (in units of the Alfven crossing time). Transition 
to non-axisymmetric equilibrium occurs at around t/r\ = 120. 

3.2.1 Analytic stability conditions 

The m ^ 2 modes always set in after the m = 1 mode and need 
not be considered here. Taking the output the fiducial simulation 
(where E/U — 1/400), described in the previous section, I apply 
Tayler's six conditions ((6) for both m = and m = 1 modes to the 
toroidal component of the field, at the three points in time shown 
in fig. [4] Looking at where in the star the conditions are met, the 
first result is that the first two conditions, a > and c > for 
both the m — and m = 1 modes are satisfied everywhere in 
the star at all points in time. In contrast, the two b 2 < Aac condi- 
tions are not met everywhere; the regions in the meridional plane 
where they are and are not satisfied this are plotted in fig. [6] for 
the case where r n = 0.47-R. Of course, this is the conclusion that 
Tayler arrived at - that a purely toroidal field is always unstable in 
at least some part of the star. The field which forms in the simu- 
lations is stable against the Tayler instability only because of the 
presence of the poloidal component. However, it is possible that 
the toroidal component can be stabilised with a weaker poloidal 
component than is actually present in the fiducial simulation. We 
can use J14t to estimate what strength the poloidal field needs to 
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Figure 6. Half of the meridional plane (the other half being identical) with 
Tn/R = 0.47 value. The thick black line is the surface of the star and the 
thin black lines are two selected poloidal field lines, showing the locations 
of the 'closed' poloidal region (where the toroidal field resides) and of the 
neutral line. The regions stable against the m = and m = 1 even in the 
absence of a poloidal field are illustrated by the green and blue lines; the 
ticks point into the region of instability. The thick pink line shows which 
area is stabilised by the radial component of the field B r . Dotted and dashed 
purple lines are contours of minimum and maximum unstable wavelengths, 
in units of the grid spacing dx fa 0.03iJ, as calculated by solving for n in 
)13t and taking A = 2n/n. 



be at any particular point, so an estimate for the minimum overall 
strength of the poloidal component is that at which this condition is 
just satisfied at every point where Tayler's criteria are not satisfied. 
It should be stressed however that this approach is only approxi- 
mate, since we have ignored curvature effects, and also because lh 
is undetermined - we shall simply assume here that ~ R. The 
area stabilised according to j I4t is also plotted in fig. [6] It is clear 
from the figure that the entire region unstable to m = 1 is stabilised 
by the poloidal field, as well as almost the entire region otherwise 
unstable to m — 0, but that two small regions are apparently still 
m — unstable. However, looking at the minimum unstable wave- 
length given by l |13t , we see that both of these areas, there is barely 
space to fit one wavelength as A m i n ~ 2dx w 0.06 J?. This explains 
the stability in the simulations at this poloidal/toroidal ratio. Note 
that the reason for expressing the minimum and maximum unsta- 
ble wavelengths in terms of the grid spacing is that we can easily 
check if the simulation is running at sufficiently high resolution - 
clearly the wavelength needs to be at least a few grid spacings to be 
properly resolved. 

If the relative strengths of the toroidal and poloidal compo- 
nents are now changed (while keeping the total energy E con- 
stant), the position of the B r stabilisation line changeQ Fig. [7J 



4 In principle the m = and m = 1 stability lines also move, but in 
the case of strong toroidal field and constant total energy, the toroidal field 
strength changes only by a small factor and the lines move very little. 
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Figure 8. Stability lines for r n = 0A7R and E p /E = 0.001 in the weak- 
field case (E/U = 1/4000). The resemblance with the E p /E = 0.01, 
E/U = 1/400 plot in fig.|7]is obvious, the only significant difference is 
that the unstable wavelengths are a factor \/T0 shorter in the weaker field 
case. See fig.|6]for an explanation of the various lines plotted. 

shows the positions of the stability lines at various ratios E p /E 
of poloidal/total energy at the three different values of r n . 

Looking first at the top row of fig. [7J(r n /72 = 0.33, E p /E = 
0.139 in the fiducial simulation), it can be seen that the E p /E — 
0.01 should be unstable, as there is an m = 1 unstable region 
on the left where the minimum unstable wavelength is only a few 
grid spacings, small enough to fit into the space available (but large 
enough to be resolved numerically in the simulations described 
below). However, the minimum wavelength looks too high in the 
m — unstable zone. The 0.018 case is more marginal but there 
may still be space to fit an unstable wavelength of 5 grid spacings 
( 0.15-R,) into the m — 1 unstable space. The 0.032 case looks also 
marginal, but 0.056 does look stable. Therefore, the critical value 
of E p /E is perhaps around 0.018 or 0.032. Repeating this analy- 
sis on the middle row of fig.|7J(r n /i? = 0.47, E p /E = 0.046 in 
the fiducial simulation) we get a stable E p /E minimum at around 
E p /E = 0.0056 to 0.01. Finally on the bottom row of fig. UJ 
(r n /R = 0.58, E p /E — 0.066 in the fiducial simulation) we see 
that practically the entire region of interest is m — unstable, 
which is stabilised by the poloidal field at around E p /E — 0.056, 
but that the m = 1 mode should be stabilised at a somewhat lower 
poloidal energy ratio, perhaps at roughly 0.01. 

These lower limits on E p /E are, as discussed in sec tion [2". 1.21 
dependent on the absolute field strength which I parametrize as 
E/U, the ratio of magnetic to thermal energies. Changing the mag- 
netic field energy in the fiducial simulation by a factor of 10 so that 
E/U — 4000 and repeating the exercise above results in stability 
at a value of E p /E a factor of 10 lower, as can be seen in fig. [8] 
where the ratios E p /E and E/U are both a factor of 10 lower than 
in the second plot of the second row of fig. [7] but the result is near- 
identical stability lines. The only significant difference between the 
two is that the unstable wavelengths are a factor \/T0 lower in the 
weaker field case. 
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Figure 7. Stability lines for various values of E p /E at three points in the fiducual simulation with E/ U = 1/400. Top row: r n = 0.33iJ and E p /E = 0.01, 
0.018,0.032 and 0.056; middle row: r n = 0.47i? and E p /E = 0.0056, 0.01, 0.018, and 0.032; bottom row: r n = 0.58i? and E p /E = 0.01, 0.018, 
0.032 and 0.056. See fig.[6]for an explanation of the various lines plotted. 



3.2.2 Simulations with different E p / E 

To verify this assessment above, it is possible to use these altered 
E p /E configurations as the initial conditions for simulations. I use 
here simulations with E/U = 1/400. First, the star is allowed 
to relax to equilibrium by keeping the magnetic field constant and 
letting the pressure field adjust until it balances the Lorentz force. 
After ten or so sound crossing times, this relaxation has taken place. 
A small white-noise perturbation is then added to the density field, 
and the magnetic field is allowed to change, i.e. the induction equa- 
tion is switched on (this defines t — in the following plots). In 
fig-IU the amplitudes of the azimuthal modes m = to 4 (an r.m.s. 
integration over the meridional plane of velocity component vg) 
are plotted against time for simulations with the following values 
of E p /E: 0.01, 0.018, 0.032 and 0.056, the initial conditions hav- 
ing been taken from the fiducial run at r n /R — 0.33 before having 
the poloidal energy ratios changed by hand. It can be seen in the 
figure that the threshold is at Ep/E w 0.032. 

This exercise was then repeated with the equilibria of different 
values of r n /R'- 0.47 and 0.58: see figs. 1101 and 1111 The former is 
easy enough to interpret: the threshold is around E p /E = 0.01. 
The small contamination we see in m = 4 from the geometry of the 
computational box is visible here but much greater in the r n /R = 



0.58 simulations, presumably because at higher r n the field has a 
stronger interaction with the sides of the computational box. Here, 
the ratio 0.056 simulation, which we know to be stable or almost 
so because the fiducial run has a ratio 0.066, has significant m = 4 
contamination but we can see that the m = 1 mode is stable at 
least for the first few Alfven crossing times, so it looks like the 
critical ratio is somewhere between 0.032 and 0.056. Surprisingly, 
the m = mode seems stable at all ratios presented in the figure. 
This is at odds to the conclusions of section[T2j] where we found 
that whilst m — 1 is stable above 0.01, the m = mode should be 
unstable up to E p / E — 0.056 or thereabouts. The reason for this 
discrepancy is not immediately clear, but may have something to 
do with effects ignored in this study such as curvature effects - the 
unstable region at this value of r n is near the equator where these 
could be important. 

At this juncture it is necessary to look at the minimum and 
maximum wavelengths plotted in fig. [7] to check that the simula- 
tions have sufficient resolution, as it is conceivable that a E p /E 
ratio which is in reality unstable could appear to be stable in a sim- 
ulation simply because the unstable wavelengths are too low to be 
resolved. Looking at the top row, where the Ep/E — 0.032 sim- 
ulation is the first which appears stable, we see that the unstable 
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Figure 9. Log amplitudes of modes ^ in ^ 4 (solid, dotted, dashed, dot- 
dashed and dot-dot-dot-dashed respectively) in simulations with r n /R = 
0.33 and E p /E = 0.01, 0.018, 0.032 and 0.056. The stability threshold 
appears to be around 0.032, below which the m = 1 mode is the dominant 
mode. 





Figure 11. Log amplitudes of modes ^ m ^ 4 io simulations with 
r n /R = 0.58 and E p /E = 0.0032, 0.0056, 0.01, 0.018, 0.032, and 
0.056. The m = mode appears to be stable at all ratios considered here, 
whilst the m = 1 mode appears to stabilise at around E p /E = 0.056 (see 
text). 



from fig.[8]that the unstable wavelengths will be badly resolved or 
not resolved at all, and therefore that significantly higher numeri- 
cal resolution would be required. At higher field strengths we are 
no longer in the physically interesting weak-field regime, as the 
structure of the star becomes significantly non-spherical. Checking 
this result numerically at the higher resolution required will be left 
for the future. It is not inconceivable of course that some other in- 
stability not considered here could become relevant at lower field 
strengths. 



Figure 10. Log amplitudes of modes ^ m ^ 4 in simulations with 
r n /R = 0.47 and E p /E = 0.0056, 0.01, 0.018 and 0.032. The stability 
threshold appears to be around 0.01 or 0.018. 

region near the axis of symmetry has a maximum unstable wave- 
length of more than 10 grid spacings (Wdx « 0.37?, about the size 
of the unstable region) and a minimum wavelength of great er than 
5dx. Now, in previous simulations (e.g. Braithwaite 2006a) it was 
found that this high-order code can resolve modes of wavelength 
8dx almost perfectly, so it seems unlikely that a A ~ 5dx mode is 
entirely suppressed. Looking at the middle and bottom rows, where 
the E p /E = 0.01 and 0.032 simulations respectively are the first 
which appear stable in the simulations, the minimum wavelengths 
in the m = 1 unstable regions are also around 5dx in both cases. 

Ideally one would now perform simulations with different 
field strengths, i.e. different values of E/U from the value 1/400 
used above, in order to check the analytic prediction that the criti- 
cal poloidal energy fraction {E p /E) CT it oc E/U. However, we see 



3.3 Stability of a predominantly poloidal field 

We have seen that an axisymmetric configuration with much 
stronger a toroidal component than poloidal can be stable. How- 
ever, we have reason to believe that the opposite is not true; we 
expect that a field with a much stronger poloidal component will 
make a transition to non-axisymmetric equilibrium. 

We can now apply the stability conditions l l 1 6t and dl7t to the 
fiducial simulation at the three points in time used in the previous 
section, i.e. when r n /R = 0.33, 0.47 and 0.58. We find that the 
thresholds in E p /E are 0.71, 0.48 and 0.69 using {[6j. How to cal- 
culate Wright's (1973) threshold is not obvious; one way is to look 
at conditions near the neutral line and use soB p /s for B p (so) and 
set so to r n , and use the value of B t on the neutral line. This gives 
thresholds of 0.91, 0.80 and 0.91. Wright points out that his stabil- 
ity condition is a necessary condition, and should be interpreted as 
a lower bound on the toroidal field required for stabilisation. 

To test the strength of toroidal field required for stabilisa- 
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Figure 12. Log amplitudes of azimuthal modes $C m $C 4 (solid, dot- 
ted, dashed, dot-dashed and dot-dot-dot-dashed respectively) in simulations 
with r n = 0.33/?. The stability threshold appears to be between 0.8 and 
0.9. 

tion, simulations analogous to those presented in section [3".2.2l were 
performed, this time with high values of Ep/E. The pertubation 
added in this case differs from the white noise used in the previous 
section: a large-scale pertubation was added to azimuthal modes 
m = 2 to 8. In fies. l 1211 13l andfT4lthe amplitudes of modes m — 
to m = 4 are plotted against time for different poloidal energy 
ratios at different values of r n . 

The thresholds in all cases seem to be at Ep/E m 0.8. In the 
r n = 0.33.R simulations (fig,ll2b. at 0.8 (and at 0.7), an oscillatory 
behaviour can be seen in the m — 2 mode, indicating stability. At 
0.9, the m — 2 mode is dominant; at 0.944 both m — 2 and m = 3 
grow (the m = 4 mode is probably just numerical contamination at 
these ratios, but will become properly unstable above some thresh- 
old in E B /E). Looking at fig.HH where r„ = 0A7R, we see that 
Ep/E — 0.7 is clearly stable; at 0.8 the m — 2 appears stable 
and then unstable - this is probably because the magnetic field is 
just below the threshold and is taken over it by continuing secular 
evolution, which weakens the toroidal component faster than the 
poloidal. It is visible in the plots that at higher poloidal energy ra- 
tios that the higher azimuthal modes become unstable. Similarly 
we see in fig.[l2](r n = 0.58ii) that the stability threshold is very 
close to 0.8. Very clear here is the stability (hence oscillations) of 
m = 3 at 0.9 but unstable growth at E p /E = 0.944. Of course, 
oscillations follow unstable growth when saturation, i.e. a new non- 
axisymmetric equilibrium, is reached. 

A perhaps more elegant way to present the output from the 
simulations is to plot a map of the field on the surface of the 
star. Figs. Q3] and [16] are such maps of B r ; all are taken from the 
r n /R = 0.47 simulations. Fig.. Q3] follows the evolution in time 
of three simulations with different E p /E ratios; by the end of the 
simulations new equilibria have been reached in the unstable cases. 
These maps confirm that E p /E = 0.7 is stable (apart from a small 
m = 4 numerical contamination), that 0.8 is stable at first and then 
succumbs to the m = 2 mode, and that the 0.9 field is unstable 
mainly to the m = 2 mode and eventually finds its way into a new 
equilibrium. Fig. [16] shows the equilibria reached in simulations 
with even higher Ep/E ratios; evidently, the higher the ratio, the 
more complex the resulting equilibrium. 



Figure 13. Log amplitudes of azimuthal modes ^ m $C 4 in simulations 
with r n = 0.47/?. The stability boundary is around 0.8. 
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Figure 14. Log amplitudes of azimuthal modes ^ m ^ 4 in simulations 
with r n = 0.58/?. 



It seems then that a magnetic field becomes unstable as the 
poloidal energy fraction exceeds about 80%, and that this thresh- 
old depends only weakly on the central concentration of the field 
(i.e. on r n ). Just above this threshold, the only unstable mode is 
m — 2, which is reassuringly well resolved by simulations at this 
resolution. At higher values of Ep/E higher modes become un- 
stable, which is exactly what we expect because higher azimuthal 
modes have to do more work against the toroidal field, in proportion 
to the value of m. Therefore, the m — 3 mode becomes unstable 
with a toroidal field 2/3 the strength of that at the m = 2 stability 
threshold, so that if the m = 2 threshold is E p /E w 0.8, then the 
m = 3 threshold will be at E p /E w 1 - (1 - 0.8)(2/3) 2 = 0.911 
(ignoring the small change in poloidal field strength resulting from 
the change in E v /E at constant E). This can be seen for instance 
in fig.[l4]and confirms the m dependence of in l !16t . 

The stability thresholds found from these simulations agree 
only qualitatively with the predictions d 1 6b and l !17t . This probably 
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Figure 15. Behaviour of the magnetic field when the poloidal component is stronger than the toroidal, near the maximum poloidal/toroidal ratio. Maps of 
the radial magnetic field B r on the stellar surface at three points in time (t/TA = 2.9, 5.8 and 10.5; plates on the left, middle and right respectively) of the 
simulations with r n /R = 0.47 and E p /E = 0.7, 0.8 and 0.9 (top, middle and bottom rows respectively). White and black represent the strongest positive 
and negative B r , and the black line shows the position of B r = 0, which is initially at the equator. The simulation on the top row shows stability; on the 
middle row we see initial stability followed by an m = 2 mode, and the simulation with E p /E = 0.9 shows a strong growth of the m = 2 mode which then 
saturates after a few Alfven crossing times as the field settles into a new equilibrium. 
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Figure 16. As fig. 1 151 but all three B r maps are at one point in time (</ta = 10.5) in three different simulations with E p /E = 0.968, 0.99 and 1.0. 
Evidently, at higher initial E p / E ratios, the non-axisymmetric equilibria reached are more complex. 



has to do with effects not included in the analysis, such as the the 
toroidal geometry of the flux tube, the fact that the poloidal field 
lines are not circular in the neighbourhood of the neutral line - in 
fact they are decidedly elliptical, and that the toroidal field falls off 
away from the neutral line. 

It is informative to think about this evolution into a new equi- 
librium in terms of the magnetic helicity, defined as J B ■ AdV 
where A is the vector potential. As the field evolves on a dynam- 
ical timescale, helicity is conserved and the field can be thought 
of as evolving into the lowest energy state for that value of helic- 
ity. Essentially, helicity can be thought of as the product of toroidal 
and poloidal fluxes, so that as we go to higher E p /E ratios, the 
helicity falls, and below some threshold the lowest energy state is 
non-axisymmetric. The initial equilibrium is essentially a twisted 
flux tube lying in a circle around the equator of the star, and the 
transition to non-axisymmetric equilibrium is a matter simply of 
stretching this flux tube into a more complex arrangement. In the 
process of stretching the tube, the toroidal component (i.e. the com- 



ponent parallel to the axis of the tube, the neutral line) is amplified 
(since the tube becomes narrower), and the poloidal component be- 
comes weaker, eventually bringing the two components to roughly 
equal strengths, because the energy minimum for a given helicity, 
i.e. for a given product of toroidal and poloidal field strengths, will 
have the two components roughly equal to each other. At higher 
E p /E ratios therefore, more stretching is required to make the two 
components equal. This can be clearly seen in fig.[T6] Now, in the 
case where E p /E = 1 the field has zero magnetic helicity and no 
amount of flux-tube stretching can result in an equilibrium. How- 
ever, there are diffusive processes at work which can either create 
helicity or split the one original flux tube into two or more tubes 
which can have helicity of different signs and which add up to zero, 
although it is likely that a lot of time will pass before any equilib- 
rium is reached and the energy of the equilibrium will be very much 
lower than the original energy. 
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4 CONCLUSIONS AND DISCUSSION 

In this paper, I have looked at the lower and upper limits on the 
fractions of energy in the poloidal and toroidal components of an 
axisymmetric magnetic field. To find these limits, I took the output 
from a simulation where a 'turbulent' initial magnetic field evolves 
into an axisymmetric equilibrium, changed the relative strengths 
of the poloidal and toroidal components by hand, and used that as 
the initial conditions for new simulations. This was supplemented 
with more analytic methods including the necessary and sufficient 
stability conditions found bv lTavlerl J1973I) (it is incidentally found 
that four of his six conditions are always met at every point in the 
star). The two methods are in broad agreement. 

The result of this investigation is that while the upper limit on 
the poloidal energy fraction E p /E is around 80%, the lower limit 
depends on factors such as the radius of the neutral line r n and 
can be between 1% and roughly 5% for a star constructed from a 
polytrope of index n — 3 (which approximates to an upper-main- 
sequence star) and where the ratio of magnetic to thermal energies 
E/U — 1/400 . This lower limit is expected to be proportional to 
the ratio E/U, so that (E p /E) clit ~ WE/U. These limits will 
also depend on other factors not explicitly explored here, such as 
the equation of state and density profile of the star, but these should 
not affect the results in more than a modest quantitative manner; 
however in a NS we might expect a lower limit of (E p /E) C rit ~ 
10 3 E/U. The upper limit found here broadly confirms what was 
expected from the analysis in Paper II and from the analyses of 
IWrightl Jl973l) and iMarkev & Tavlej Jl974l) . who found that the 
toroidal field must be at least about a quarter of the strength of the 
poloidal field. The lower limit on E p /E had not been examined 
previously. 

The question of what ratios we actually expect to find in nature 
has not yet been answered. It will depend on the state of the mag- 
netic field left over from the convective protostellar phase (whether 
we are looking at main-sequence stars, white dwarfs or neutron 
stars) and on the subsequent secular evolution. In the fiducial sim- 
ulation described, E p /E does fall to around 0.046 and is below 
0.1 for most of the period of diffusive evolution. Given that a star 
is generally strongly differentially rotating when it is formed, and 
that any seed field will be wound up and predominantly toroidal, 
it seems plausible that the eventual equilibrium could have rather 
low Ep/E. A proper study of the effect of initial conditions on the 
resulting equilibrium will be left for the future. 

It remains to looks at some of the implications of these re- 
sults. The magnetic field of a star can have various effects on its 
appearance and behaviour. In upper main-sequence stars and white 
dwarfs, we can directly observe the field on the surface via the Zee- 
man effect, and it turns out that many, if not most, do have roughly 
axisymmetric configurations. In neutron stars, we measure the spin- 
down and infer from that the dipole component on the surface. Be- 
low the surface of the star, there could be a deeply buried field (with 
low r n and very low $ sur f /$ p ) and/or a very strong toroidal com- 
ponent, both of which could in effect 'hide' magnetic energy from 
view. This is of obvious interest in the study of magnetars, neutron 
stars with observed dipole fields of 10 14 ~ 15 gauss whose emission 
in X rays and 7 rays is powered by the decay of the magnetic field. 
The possibility that these stars could contain a large quantity of 
energy hidden from view could explain the large energy output of 
these objects. A field of 10 15 G contains around 2 10 47 erg in mag- 
netic energy, and flare has been observed which emitted a tenth of 
that quantity in less than a second. According to the standard flare 
model, some slow evolution of the magnetic field in the core re- 



sults in stress build-up in the crust, which eventually results in the 
crust cracking and a release of energy, but it is difficult to imagine 
this mechanism releasing a large fraction of the magnetic energy 
during any one event. A field which is stronger in the core of the 
star than the poloidal component we see on the surface is one so- 
lution to this problem. Another possible solution is that the field is 
very non-axisymmetric, and that the dipole measured is an under- 
estimate of the average field strength on the surface. Any of these 
are also potential solutions to the phenomenon of neutron stars with 
very different observational properties but which occupy the same 
region on the P — P diagram and have therefore the same dipole 
field strength. 

Another effect that the magnetic field has is to deform the 
star's mass distribution. In the light of results presented above, 
it seems likely that a star has a predominantly toroidal field, de- 
forming the star into a prolate shape. This causes the star to un- 
dergo torque-free precession and the damping of this precession 
will cause the star's magnetic axis to tend towards orthogonality 
with the rotation axis. In neutron stars this could potentially happen 
fairly quickly, faster than the orientation of the axes can be changed 
by spindown torque, which acts on the spindown timescale. In the 
millisecond magnetar models, the proto-NS rotates at some signif- 
icant fraction of break up and this results in a powerful dynamo 
which creates a strong magnetic field during the first hundred sec- 
onds. Since it seems very likely that the equilibrium formed after 
this dynamo switches off is predominantly toroidal, the star should 
flip over and emit gravitational waves observable with the next gen- 
eration of detectors at least as far away as the Virgo cluster. How- 
ever, the spindown timescale of a NS rotating with P — 1ms and 
B — 10 15 G is only an hour, so the flipping-over mechanism has 
to work fairly fast. Since the free precession period of such a star 
would be P/e = 1000s, where e ~ 10 -6 is the ellipticity induced 
by the magnetic field of this magnitude, the damping of the pre- 
cession would have to be very efficient indeed unless the magnetic 
equilibrium could be produced already at an angle to the rotation 
axis. More work is required to address these issues. 
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